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The Social System Identification Problem 

Hoi-To Wai^^, Anna Scaglionel^, Amir Leshem^ 


Abstract —The focus of this paper is modeling what we call 
a Social Radar, i.e. a method to estimate the relative influence 
between social agents, by sampling their opinions and as they 
evolve, after Injecting in the network stubborn agents. The 
stubborn agents opinion is not influenced by the peers they seek 
to sway, and their opinion bias is the known input to the social 
network system. The novelty is in the model presented to probe 
a social network and the sointlon of the associated regression 
problem. The model allows to map the observed opinion onto 
system equations that can be nsed to infer the social graph and 
the amonnt of trnst that characterizes the links. 

1. Introduction 

Recently, the rapid growth of online social media such 
as Facebook, Twitter has generated a lot of interests in 
researches on social networks. Importantly, it has provided 
a platform for researchers from multiple disciplines, rang¬ 
ing from social science to statistical physics, to study and 
understand the behavior of the human society. In the con¬ 
trols area, due to the close ties with decentralized robots 
coordination problems, several works focuses on modeling 
the opinion dynamics/exchanges (see e.g. [l]-[8] and refer¬ 
ences therein). Arguably, understanding the opinion dynam¬ 
ics through which individuals seek to share information and 
agree is one of the most important social studies. In light of 
this, extensions to the basic opinion dynamics model are also 
popular, e.g., [9]-[14]. Also, in the recent developments, the 
notion of controllability and observability has been extended 
to the context of complex networks [15]-[18]. 

The common definition of a social network focuses on the 
social graph component [19], where we represent individuals 
(social agents) as nodes and the friendships between them 
as edges. Knowing the social graph alone is insufficient for 
understanding social networks. In particular, individuals may 
exhibit different degree of trust in their neighbors. There is 
strong trust among the close friends and weaker trust between 
individuals without mutual interests [20], [21]. Identifying 
the social system allows us to predict the behavior of 
individuals in a social network in times of decision making. 
The difficulty is that, while the interaction between agents 
may be evident, the trust between them and the impact an 
interaction has on another agent is not directly observable. 

The focus of this study is identifying what we call a social 
system, which encompasses both the social graph and the 
set of trusts between individuals. We define the state as 
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the instantaneous opinions; and the output as the observed 
opinions after a certain period of time of discussions. The 
social system is modeled as an endogenous system of opinion 
updates that follows exogenous stimuli that we observe. 

We are interested in solving the social system identification 
(SSI) problem, which could be viewed as the development 
of a Social Radar since the idea is injecting test signals 
into the system and observing its output [22], much like 
in a traditional Radar. An inverse problem is then solved 
by gathering observations that are tied to the input opinions 
and output opinion pairs. The challenge is that the latter are 
endogenous. In the literature, a related issue is the inference 
problem of graphical models, e.g., [23]-[25]. The latter 
assumes that the effect of the social system is manifested 
in the correlation of the neighboring opinions, instead of 
the opinion dynamics. Indeed, our model is similar to that 
considered in [26]-[28]. Particularly, the methods proposed 
in [27], [28] consider the case with non-linear dynamics. 
However, the methods require knowing precisely when the 
opinion of an agent has impacted that of another, which we 
claim is unrealistic, given that what happens in people brains 
is not visible and tracking their impact on an individual 
would require testing the opinion of the neighborhood again, 
which is an unnatural way of communicating. 

In contrast, our work assumes only partial knowledge on 
the social graph and that opinion updates have occurred at 
unknown times due to unknown stimuli. All that we observe 
are noisy versions of the agents opinions. This motivates us 
to use steady state models as an approximation for what ties 
the opinions we sample over time. However, to identify the 
system we need, therefore, to prevent trivial consensus. Our 
idea is to introduce a set of stubborn agents, i.e., agents who 
are not swayed by other opinions [11], [12], [29], [30], into 
the social network. The stubborn agents serve as ‘probes’ 
inserted on the social network that injects input to a social 
system. Indeed, the stubborn agents change the terminal 
behavior of the opinion dynamics and reveal the social 
system in the form of an underdetermined linear system. 

The contributions of this paper are two-fold. Firstly, we 
formulate the SSI problem under the presence of stubborn 
agents, and provide a set of conditions for identifiability. 
Secondly, we consider the random opinion dynamics model 
and propose an estimator for the (ensemble) mean of terminal 
opinions. The mean square convergence of such estimator is 
proven. We provide numerical results to verify our findings. 

Notations: The Kronecker product is 0 and vec(-) as the 
vectorization operator. Moreover, Diag : M" —> and 

diag : > K" are defined as the diagonal operators on 

square matrices and vectors, respectively. 



II. System Model 

The social network has an associated graph G = (V, E) 
where the vertexes are the set of agents V = \ri\ = 
{ 1 , 2 ,n} and social system we want to explore is the tuple 
S = {E, W) where W is the row-stochastic matrix of the 
trust coefficients. We assume: 

Assumption 2.1: The trust matrix W satisfies Wij > 0 
and/or Wji > 0 if and only if ij G E. Wu is the self-trust. 
Our goal is to identify W by observing opinions whose 
dynamics are consistent with S. Specifically, we assume 
that the agents are shaping opinions over a certain issue 
s S N. Initially, each agent holds an opinion (belief) on 
a discrete random variable 0 s S i.e. the 

p.m.f. Xi( 0 ;s) = (p( 6 »i|si,s),p( 6 ' 2 |si,s), 
where p{9j\si^s) is the ith agent’s belief on the event 0 s = 9j 
and Si s is the private information agent i has before its 
interactions. Notice that p{0m\si^s) = 1 — l^Xi(0;s). The 
beliefs are forged by the DeGroot’s model [20]: 

Xi{t + l;s)= (1) 

j&Mi 

where Afi is the set of neighbors of i and Wij (f) is the (/, j)th 
element of the row stochastic matrix W{t). We assume that 
W (t) is i.i.d. and drawn from a p.d.f. satisfying E{VF (f)} = 
W. The opinion dynamics can be described as 

x{t + l;s) = W{t)x{t;s), (2) 

where x{t\s) = (xi(i; s),..., x„(f; s))^ stacks the 1 x m 
vectors to form an n x m matrix. Notice that the above 
dynamics includes the randomized models in [31], [32] as 
special cases. Moreover, what is active at a given time is 
random and the trust matrix W is embedded in the opinion 
dynamics Our observations are actions/ratings that an 
agent performs and that the social network is exposed to 
(e.g., ‘liking’ a post on Facebook). We assume it is possible 
to take a noisy snapshot of the opinions at time tp. 

y{ti-,s) = x{tp,s) + n{ti-s), (3) 

where nfps) contains i.i.d. noise samples with bounded 
variance It is due to the fact that we do not have direct 
access to the opinion or are using incomplete or outdated 
information about it. Eq. (|^ & ([^ give a linear system 
representation for the social system S with the state being 
the opinions. 

We define the social system identification (SSI) problem 
as the task of inferring W from a set of measurements 
y{Ts;s) = {y{U]s)}t ^gT; (and over a number of issues 
s = The general SSI problem is challenging to 

solve for several reasons. For example, we see that W is 
hidden in the random model Q and Q; also, from © it is 
not even possible to retrieve W{t) from the samples since 
m n, i.e., the linear equation is rank-deficient. We see that 
additional prior knowledge must be incorporated to develop 
a tractable SSI method. 

There are a few prior studies on the SSI problem. Most 
closely related to ours is the work in [26]-[28]. In particular. 



Fig. 1. Using stubborn agents as ‘probes’ on social systems 

[26] considers the same model as ours. The authors assume 
that the set Ts is consecutive, e.g., Ts = {ti,ti + 
and the trast matrix is static W (t) = W with known sparsity 
pattern. In [27], [28], the authors consider a nonlinear dy¬ 
namical system and applied compressed sensing to infer the 
network topology from samples marked with time stamps. 

The assumptions made in [26]-[28] may be restrictive for 
social networks as the rate of interaction is unknown for 
the latter, therefore the time stamp information cannot be 
obtained accurately. Our idea is to introduce a set of stubborn 
agents as probes', see Fig. The developed SSI method 
requires only partial knowledge on the topology E and is 
applicable to the scenario with time-varying trust matrix. 


III. Stubborn Agents 


In social networks, stubborn agents are those members 
who place zero trust on their neighbors. In our framework 
we assume there are rig < n stubborn agents in the social 
network that we know of or control. We index them by 14 = 
[ris]. The resulting trust matrix W{t) is: 


W{t) = 


In, 0 \ 

Bit) Bit) ) 


(4) 


and similar structure can be found in W. We assume: 


Assumption 3.1: The support of B, £b = {ij ■ Bij > 
0} = E{V,Vs), is known. Moreover, each agent in V has 
non-zero trust on at least one agent in 14 . 

The assumption asserts that each non-stubborn agents is 
influenced by at least one stubborn agents. As the graph 
G is connected. Assumption m] implies that the principal 
submatrix D satisfies \\D lb < 1- 

We demonstrate the specific structure in 0 gives rise to 
a set of equations that allows a tractable solution to SSI. To 
begin with, we assume that the opinion exchange is static, 
i.e., W (t) = W for all f; this assumption will be relaxed 
later in Section |IV] Observe the following: 

Observation 3.2: [11] Under the assumption that 

W(t) = W. Consider (|^ by setting t —)■ oo, we get: 

lim x{t', s) = W°°x{0', s), (5) 


where 




Ins _ 0\ 

iI-D)-^B 0 J ■ 


(6) 


Consequently, non-stubborn agents’ opinions satisfy: 


lim (I 

t—¥00 


-^)^[n]\[ns] (^5 '5) 


lim Bxinpit',s) 

t—¥oo 


( 7 ) 






where we have defined a;[„^](t;s) = (xi(t; s))ig[„j. In fact, 
the right hand side of (|7]i can be replaced by Sa:[„^](r;s) 
for any r G N as the stubborn agents never change opinions. 

The final opinions are driven by the initial opinions at the 
stubborn agents. Importantly, Eq. 0 gives a set of linear 
equations that characterizes the social system S. To this end, 
we can define an estimator for limt_>.oo x{t; s): 


where the third equality is due to SI + ofF(£))l = 1 — 
diag(r>). Q.E.D. 

We define an equivalent relation ^ as: 


{B, D) ~ (S', D') : 3A > 0 s.t. ([T^ holds. 


(15) 


{'Ts\s) = 


U^Ts 


where the sampling set Ts is defined as 

Ts = {ti-.i& [| 7 ^|], U > To, V i}, 


( 8 ) 


(9) 


where To ^ 0. Notice that the time indices ti nor their orders 
are not required in the computation of (|^, i.e., we do not 
need to know the exact time in which opinion updates have 


The relative trust weights, defined as = Di^jiX — Da) 
and Blj = BijjiX — Du), is preserved for all tuples be¬ 
longing to the same equivalence class. Moreover, there exists 
(S', S') with diag(S') = 0 such that (S, S) - (S', S'|] 
Hence, we adopt a pragmatic approach to remedy the scaling 
ambiguity by fixing diag(S) = 0. 

The next ingredient is the fact that a social graph has 
typically a sparse set of links between agents, i.e., S 
is sparse. This motivates us to consider the following £q 
minimization [33] problem: 


occurred when sampling y{ti;s). 

In the case of static exchange, it is easy to check that 

min 

B,D 

I|off(r>)||o 

(16a) 

the estimator ([ 8 ]) is consistent as \Ts\ —t 00 and To —)■ 00 ; 

S.t. 

\\{I-D)Y -BZ\\l<e, 

(16b) 

see Section IV for a further discussion of its convergence 
properties. Notice that unlike [26], we do not require the 


(B D)l = 1, > 0, B>0, 

(16c) 

sampling set Ts to be composed of consecutive indices. 


Bij > 0, V ij G £b, Du =0, V 

(16d) 


Consider collecting the estimate 0 for K issues (i.e., s = 
1,..., K) into data matrices, we have the linear equation: 

{I -D)Y = BZ + N, 


( 10 ) 


where 


(n—ns ) X Km 


( 11 ) 


denote the data matrices for the opinions at the normal agents 
and stubborn agents, respectively, and N is the additive noise 
with variance cr^ that captures the estimation error from ([^. 

A. Identifying the Social System 

Our next endeavor is to formulate the respective inverse 
problem for SSI. In particular, our goal is to find the tuple 
{B,D) that satisfies the system of equations: 

[I -D)Y ^BZ T N, (BD)1 = 1, B,D>0. (12) 

We observe the following: 

Lemma 3.3: Assume AT = 0, there exists a tuple {B,D) 
that satisfies ( |T2l i. The tuple [B',D') also satisfies ( fl^ with 

B' = AB, off(D') = Aofr(D), (13a) 

diag(£)') = 1 - A(B1 + ofr(£))l), (13b) 

where off(D) denotes the square matrix with only off- 
diagonal elements in D and A is any non-negative diagonal 
matrix such that diag(D') > 0. 

Proof: The existence of [B, D) is ensured by picking 
B = B,D ~ D. It is also obvious that the second equation 
in ( [T 2 I 1 is satisfied by {B', D') for an arbitrary diagonal 
matrix A. For the first equation in ( [T^ , we observe that 

D'Y + B'Z = (Diag(diag(I9')) + oS{D'))Y -f B'Z 
= Y- A((Diag(Bl -f ofr(D)l) - ofr(r>))y - BZ) 

= Y- A((Diag(l - diag(i9)) - ofr(r>))y' - BZ) 

= Y-A{Y-DY-BZ)=Y, 


where e > 0 is a regularization parameter that depends on 
an- Notice that the last constraint is due to Assumption 2.1 
and the prior knowledge on £b- 

Problem •EH is non-convex. However, Problem GD can 
be readily convexified by replacing the Iq norm in the 
objective function by an norm: 


min ||off(£))||i s.t. Eq. ( |16b| ) to (|16d|. (17) 

B.n 

The convexified problem can be solved using off-the-shelf 
softwares, e.g., CVX [34]. 

Next, we study conditions under which ( [T 6 | ) can identify 
the social system S. Intuitively, we see that the identifiability 
condition depends on the number of stubborn agents and 
the sparsity of the trust matrix D. Importantly, in the case 
with optimized placement of stubborn agent]^ we have 
the following condition, whose proof can be found in the 
extended version of this paper [35]. 

Theorem 3.4: Let n > Ug and the support of S G 
be constructed such that each row has d non-zero ele¬ 
ments, selected randomly and independently. Define bmin = 
Tj^j, bmax — “ /dci, 

(3' = P — d/n and 5 = 1 — 1/(5 — 1). If 


H{a)+p'H{a/P') - 
a\og{P'/a) 

(2d - 3) - 1 - 2b„,ax > 0, 


d > max jd, 1 -f 


(18) 

(19) 


where H{x) is the binary entropy function, and ||<i’'||o < 
an/2 for all i, where w/ is the ith row of D’", then as 
n — )■ 00 , solving ( [T^ yields {B*,D*) — {B^,D^). 

*The con'esponding A can be found as A = Diag(l/di,..., l/dn—n^) 
with d = B1 + ofF(£))l. 

^Such is possible in a controlled experiment setting, where the stubborn 
agents can be controlled to influenced a sub-group of ordinary agents. 







Using ( [T 8 ] i it is possible to derive a lower bound f3{d, a) 
on P that depends on d, a and thus the number of stubborn 
agents required. In fact, P{d, a) is a decreasing function in 
d. We note from the proof of the theorem that there is a 
tradeoff between d and the probability of successful recovery. 
As such, the parameter d has to be chosen judiciously. In 
addition, condition ( [T^ requires the stubborn agents to be 
sufficiently influential to the non-stubbom agents. 


IV. Randomized Models with Stubborn Agents 

This section considers a general model of with ran¬ 
domized opinion exchange where W (t) is time varying and 
i.i.d. with mean W. Under this setting, the limit equation 
(0 only holds in expectation, i.e., 

lim (I - s)|a;(0;s)} 

— ( 20 ) 

= lim s)|a;( 0 ; s)}. 


Notice that the expectation is taken over the ensemble of 
sample paths of {W (()}t. In practice, computing the expec¬ 
tation requires the social network to ‘repeat’ the discussion 
on the same issue. Obtaining lim(_>.oo E{a:(f; s)|a:(0; s)} can 
be difficult in terms of implementation. 

Interestingly, it has been observed that in a randomized 
opinion exchange model, the introduction of stubborn agents 
leads to a behavior known as opinion fluctuation [29], [30]. 

Observation 4.1: If ris > 2 and the opinion exchange 

model is random, then x{t -b 1 ; s) x{t] s) almost surely. 

Our goal is to derive an estimator for 

limt_>oo E{a;(f; s)|a;(0; s)} that relies on the samples 
from a single issue s only. From Observation 4.1 a natural 
design is to consider an estimator that averages over the 
temporal samples, i.e., Eq. ([ 8 ]l. We have [36]: 

Theorem 4.2: Consider the estimator in ([^ 
with the sampling set Ts- Denote x{oo] s) — 

limt_>,oo E{a;(f; s)|a;(0; s)} = W x{0-,s) and assume 

that ||I 9||2 = E{|j£)(t)|| 2 }. If To —)• oo, then 

1 ) the estimator ([^ is unbiased: 

E{®(7i; s)|a;( 0 ; s)} = x{oo] s). ( 21 ) 


2 ) the estimator is consistent: 

lim E{||i(7^;s)-®(oo;s)|||.|a;(0;s)} = 0. (22) 

|Ts|— s-oo 

It follows that the method in Section |III-A| can be applied. 

Note that the convergence of 0 is related to the ergodicity 
of the random process (|^. For instance, [4] has studied the 
convergence of the ergodic mean of an opinion exchange 
model with external input. To our knowledge, our result is the 
first for randomized opinion exchange with stubborn agents. 


A. Proof of Theorem \4.2\ 

We first prove that the estimator is unbiased. Consider the 
following chain: 

E{*(Ts;s)|a;(0;s)} =^ E{y((i; s)|a:(0; s)} 

' UaTs 

= E W''x{0- s) = W^x{Q- s), 

( 23 ) 


where we have used the fact that To ^ oo and f > To for 
all ti in the last equality. 

Next, we prove that the estimator is asymptotically con¬ 
sistent, i.e., ( |2^ . Without loss of generality, we let ti < 
t 2 < ... < (| 7 ;| as the sampling instances. The following 
shorthand notation will be useful: 


$(s, t) = W{t)W{t - 1)... W(s -b l)W{s), 


(24) 


where t > s and is a random matrix. Our proof 

involves the following lemma: 

Lemma 4.3: When |f — s| —> oo, the random matrix 
4>(s,() converges almost surely to the following: 


lim = 

|t —s|—>-oo 


I 0 

B{s,t) 0 


(25) 


where ... D{q))B{q) is bounded 

almost surely. 

The proof is in Appendix We consider the following: 
E{|1*(7;; s) - ®(oo; s)|||.|a:( 0 ; s)} = 

= E {y{U]s)-x{oo]s)) ^|a;( 0 ;s)|. 


tidTs 


(26) 


Recall that y(ti]s) = x(ti; s) -b n(ti; s) and the noise 
term n{ti;s) is independent of W{t) for all t. The above 
expression reduces into: 

2 


StiGTa p}- 


(27) 


It is easy to check that the latter term vanishes when |7^| — 
00 . We thus focus on the former term. 

E {x{f-,s)-x{oo]s)) ^|a;(0;s)| 




1 


ugt. 


ir.p 

1 


E 




ueTs 


(28) 


E{Tr(Ha;(0; s)x{0; s)^)}, 


where 




(29) 


t j ^'l~s 


tiGTs 


Expanding the above product yields two groups of terms 
— when ti = tj and when ti f tj. When ti — tj, using 
To ^ oo and Lemma [43] it is straightforward to show that: 


||E{($(0,f,) - - W°°)}\\ < C, (30) 

for some constant C < oo. As a matter of fact, we observe 
that the above term will not vanish at all. This is due 
to Observation 4.1 the random matrix does not 

converge in mean square sense. 

For the latter case, we assume tj > f. We have 


= + l,t,)^iO,f) - - TU'"). 


( 31 ) 
























Taking expectation of the above term gives: 

E{($(0, U) - ($(0, U) - W°°)}, 

(32) 

where we have used the fact that + l,tj) is indepen¬ 
dent of the other random variables in the expression and 
= W°° for any f > 0. Now, notice that 

= W°° + (33) 


for some 0 < A = Xmax{D) < 1. This is due to the fact 
that D is sub-stochastic. 

As To ^ oo and by invoking Lemma 4.3 the matrix 
(4>(0,L) — W ) has almost surely only non-empty entries 
in the lower left block. Through carrying out the block matrix 
multiplications and using the boundedless of $(0,L), it can 
be verified that 


||E{($(0,f,) - Ty“)^($(0,L) - ■W^“)}|| < 0(A‘^-‘-). 

(34) 

Combining these results, we can show 

E{Tr(5a;(0;s)a;(0;sn} ^ 

ir.p > 

(35) 

for some C < oo. Notice that min^ \tk+i — ffc| > i and 
the terms inside the bracket can be upper bounded by the 
geometric series Consequently, the mean 

square error goes to zero as \Ts\ oo. The estimator ([^ 
converges in the mean square sense and is thus consistent. 

Remark 4.4: From ([3^, we observe that the upper bound 
on mean square error can be minimized by maximizing 
minij \ti — tj\. When the samples y{Ts] s) are taken from 
a finite interval [Tmax] \ [To], Tmax < cxo and |7^| < oo, the 
best estimate can be obtained by using sampling instances 
that are drawn uniformly from \Tmax\ \ [To]. 


V. Numerical Results 

This section provides numerical results for the perfor¬ 
mance of SSI. Two simple scenarios on the static and 
randomized model are considered. 

The social graph G — (V, E) is generated as an Erdos- 
Renyi graph with connectivity pe = 0.15. We fix the number 
of normal agents at n — ns = 50. We assume m = 1 and the 
number of issues to be discussed is K = 100. For the static 
model, the trust matrix W is first generated with uniformly 
distributed entries, which are then normalized to satisfy row- 
stochasticity as well as the sparsity pattern according to E-, 
cf. Assumption 2.1 For the randomized model, we have 
adopted the randomized broadcast gossip exchange model in 
[32]. In particular, at each time, a random agent wakes up and 
broadcast his/her opinion to the neighbors. The neighbors 
then mix the opinion with the weight 7 = 1/2; see [32]. 

In light of Lemma 3.3 we compare the relative trust that 
an agent has on his/her neighbors. In particular, we evaluate 
the error in estimating the relative trust matrix D as 
- D,i), - Du) and 0^ = 0 

for all i. The normalized mean square error (MSE) for D is 
/{Y^^ AdT)"^) (and similarly for B). 



Fig. 2. Average performance of social system identification under the static 
model: (Left) the relative MSE in estimating D and B (Right) the support 
recovery error — \{ij : Dij 7 ^ 0 , Dij = 0 , or, Dij = 0 , Dij 7 ^ 0 }|. 



Fig. 3. Identifying a social system under the randomized gossip exchange 
model: (Left) the normalized trust matrix D of the actual social sys¬ 
tem. (Right) the normalized estimated trust matrix . The dai'kness of 
the dot represents the amount of trust between a pair of agents. 


We first evaluate the SSI performance under the static 
model with Monte-Carlo simulation. For each we av¬ 
erage over 100 instances of social graphs to evaluate the 
normalized MSE. The noise is assumed to be zero. As 
such, we can set To = 10"* and |7/| = 1 for the estimator 
We compare the normalized MSE in terms of the relative 
trust matrices, against the number of stubborn agents Ug. The 
simulation results are shown in Fig. As seen, the system 
identification performance gradually improves as ris grows. 
Moreover, the performance is significantly better when the 
subgraph between stubborn and non-stubbom agents is con¬ 


structed as a random regular bipartite graph, cf. Theorem 3.4 
Notice that the theorem’s condition requires Ug > 38. 

The next simulation example considers the case with 
random opinion exchanges. We focus on one instance of the 
social graph generated and we set the number of stubborn 
agents to Ug = 30. We conduct the test when 17^ is generated 
as a random non-regular bipartite graph of connectivity 
Pg = 0.15. For the estimator proposed in Section IV we set 
To = 10^ and |7/| = 3 x 10^, where sampling instances ti in 
Tg is uniformly drawn from {10^ -f 1,..., 10®}. The noise 
variance is cr^ = 10“^ and we set e = '/K x 1.65 x 10“^ 
in The estimated social system is depicted in Fig. 
Notice that in this case, the normalized MSE is evaluated as 
7.01 X 10-2. 

We observe that the identified social system is close to the 
actual social system. However, some links with weak trusts 
























can also be found in the estimate D^. This is possibly an 
artefact from the estimation of a;(oo; s) using 

VI. Conclusions 

We have defined the SSI problem for identifying both the 
social graph and mutual trusts between individuals in social 
networks. The system identification is achieved via the in¬ 
clusions of stubborn agents and conditions for identifiability 
are proven. We have proposed a consistent estimator for the 
ensemble mean opinions in randomized gossip model. 

This work paves a key stone towards developing a Social 
Radar that estimates the relative influence between individu¬ 
als. Future directions will include developing an efficient and 
parallelizable solution method for solving ([T^ and tightening 
the necessary condition for identifiability. 

Appendix 

A. Proof of Lemma \4.3\ 

We first establish the almost sure convergence of 
D{t)D{t — 1)... D{s) to 0. Define 

P{s,t)^\\D{t)D{t-l)...D(s)h, (36) 

and observe the following chain 

E{^(s,f)|/3(s,f- l),...,/3(s,s)} 

< E{\\Dit)\\2\\D{t - 1)... D{s)\\2ms,t- 1)} (37) 

= ^{\\Dit)\\ 2 }P{s,t- 1) < c/3{s,t- 1), 

where c = \\D lb < 1 due to Assumption 
sure convergence of follows from [37, Lemma 7]. 

Now, expanding the multiplication ( |24l l yields: 

The desired result is achieved by observing D{t)... D{s) —>■ 
0 and B{s,t) is bounded almost surely. 

References 

[1] D. Acemoglu, M. Dahleh, I. Lobel, and A. Ozdaglar, “Bayesian 
learning in social networks,” in Review of Economic Studies, vol. 78, 
2010, pp. 1201-1236. 

[2] D. Acemoglu and A. Ozdaglar, “Opinion Dynamics and Learning in 
Social Networks,” Dynamic Games and Appln., vol. 1, pp. 3-49, 2010. 

[3] P. Jia, A. Miitabatabaei, N. E. Friedkin, and F. Bullo, “Opinion 
Dynamics and the Evolution of Social Power in Influence Networks,” 
SIAM Review, pp. 1-27, 2013. 

[4] C. Ravazzi, P. Frasca, R. Tempo, and H. Ishii, “Ergodic Randomized 
Algorithms and Dynamics over Networks,” IEEE Trans. Control of 
Network Sys., pp. 1-11, to appear. 

[5] A. Tahbaz-Salehi and A. Jadbabaie, “A necessaiy and sufficient 
condition for consensus over random networks,” IEEE Trans. Autom. 
Control, vol. 53, no. 3, pp. 791-795, 2008. 

[6] B. Touri and A. Nedic, “On ergodicity, infinite flow, and consensus in 
random models,” IEEE Trans. Autom. Control, vol. 56, no. September, 
pp. 1593-1605, 2011. 

[7] R. Carli, F. Fagnani, P. Frasca, and S. Zampieri, “Gossip consensus 
algorithms via quantized communication,” Automatica, vol. 46, no. 
September, pp. 70-80, 2010. 

[8] V. D. Blondel, J. M. Hendrickx, A. Olshevsky, and J. N. Tsitsiklis, 
“Convergence in multiagent coordination, consensus, and flocking,” in 
Proc CDC-ECC '05, vol. 2005, 2005, pp. 2996-3000. 

[9] G. Deffuant, D. Neau, F. Amblard, and G. Weisbuch, “Mixing beliefs 
among interacting agents,” in Adv. Compl. Syst., vol. 3, 2000, pp. 87- 
98. 


[10] R. Hegselmann and U. Krause, “Opinion dynamics and bounded 
confidence models, analysis and simulations,” in Journal of Artificial 
Societies and Social Simulation, vol. 5, 2002. 

[11] M. E. Yildiz and A. Scaglione, “Computing along routes via gossip¬ 
ing,” IEEE Trans, on Signal Process., vol. 58, no. 6 , pp. 3313-3327, 
2010 . 

[12] E. Yildiz, D. Acemoglu, A. Ozdaglar, A. Saberi, and A. Scaglione, 
“Discrete opinion dynamics with stubborn agents,” SSRN eLibrary, 
2011 . 

[13] V. D. Blondel, J. M. Hendrickx, and J. N. Tsitsiklis, “On Krause’s 
consensus formation model with state-dependent connectivity,” IEEE 
Trans. Autom. Control, vol. 54, pp. 2586-2597, 2009. 

[14] L. Li, A. Scaglione, A. Swami, and Q. Zhao, “Consensus, polarization 
and clustering of opinions in social networks,” IEEE J. Sel. Areas 
Commun., vol. 31, no. 6 , pp. 1072-1083, 2013. 

[15] T. Wang, H. Krim, and Y. Viniotis, “Analysis and Control of Beliefs 
in Social Networks,” IEEE Trans, on Signal Process., vol. 62, no. 21, 
pp. 5552-5564, 2014. 

[16] M. Doostmohammadian and U. A. Khan, “Graph-theoretic distributed 
inference in social networks,” IEEE J. Sel. Topics Signal Process., 
vol. 8 , pp. 613-623, 2014. 

[17] Y.-Y. Liu, J.-J. Slotine, and A.-L. Barabasi, “Controllability of complex 
networks.” Nature, vol. 473, pp. 167-173, 2011. 

[18] -, “Observability of complex systems.” PNAS, vol. 110, no. 7, pp. 

2460-5, 2013. 

[19] M. O. Jackson, Social and Economic Networks. Princeton, NJ, USA: 
Princeton University Press, 2008. 

[20] M. DeGroot, “Reaching a consensus,” in Journal of American Statist- 
cal Association, vol. 69, 1974, pp. 118-121. 

[21] N. E. Friedkin and E. C. Johnsen, “Social influence and opinions,” 
The Journal of Mathematical Sociology, vol. 15, pp. 193-206, 1990. 

[22] T. Soderstrom and P. Stoica, System Identification. Upper Saddle 
River, NJ, USA: Prentice-Hall, Inc., 1988. 

[23] M. J. Wainwright and M. I. Jordan, “Graphical models, exponential 
families, and variational inference,” Foundations and Trends in Ma¬ 
chine Learning, vol. 1, no. 12, pp. 1-305, 2008. 

[24] A. Anandkumar, V. Y. F. Tan, F Huang, and A. S. Willsky, “High¬ 
dimensional structure estimation in ising models: Local separation 
criterion,” Annals of Statistics, vol. 40, no. 3, pp. 1346-1375, 2012. 

[25] G. Bresler, “Efficiently learning Ising models on arbitrary graphs,” 
p. 20, 2014. [Online]. Available: http://arxiv.org/abs/1411.6156 

[26] A. De, S. Bhattacharya, P. Bhattacharya, N. Ganguly, and 
S. Chakrabarti, “Learning a Linear Influence Model from Transient 
Opinion Dynamics,” CIKM '14, pp. 401^10, 2014. 

[27] M. Timme, “Revealing network connectivity from response dynamics,” 
Physical Review Letters, vol. 98, no. 22, pp. I^, 2007. 

[28] W.-X. Wang, Y.-C. Lai, C. Grebogi, and J. Ye, “Network Reconstruc¬ 
tion Based on Evolutionary-Game Data via Compressive Sensing,” 
Physical Review X, vol. 1, no. 2, pp. 1-7, 2011. 

[29] D. Acemoglu, G. Como, F. Fagnani, and A. Ozdaglar, “Opinion 
Fluctuations and Disagreement in Social Networks,” Mathematics of 
Operations Research, vol. 38, no. 1, pp. 1-27, Feb. 2013. 

[30] W. Ben-Ameur, P. Bianchi, and J. Jakubowicz, “Robust Average Con¬ 
sensus using Total Variation Gossip Algorithm,” in VALUETOOLS, 
2012, pp. 99-106. 

[31] S. Boyd, A. Ghosh, B. Prabhakar, and D. Shah, “Randomized gossip 
algorithms,” IEEE Trans. Inf Theory, vol. 52, no. 6 , pp. 2508-2530, 
Jun. 2006. 

[32] T. C. Aysal, M. E. Yildiz, A. D. Sarwate, and A. Scaglione, “Broadcast 
gossip algorithms for consensus,” IEEE Trans, on Signal Process., 
vol. 57, no. 7, pp. 2748-2761, 2009. 

[33] E. Candes and T. Tao, “Decoding by Linear Programming,” IEEE 
Trans. Inf. Theory, vol. 51, no. 12, pp. 4203-4215, Dec. 2005. 

[34] M. Grant and S. Boyd, “CVX: Matlab software for disciplined convex 
programming, version 2.1,” http://cvxr.com/cvx Mar. 2014. 

[35] H.-T. Wai, A. Scaglione, and A. Leshem, “Active Sensing of Social 
Networks,” submitted to IEEE Trans. Sig. and Inf. Proc. over Net¬ 
works, 2015. 

[36] S. M. Kay, Fundamentals of Statistical Signal Processing: Estimation 
Theory. Upper Saddle River, NJ, USA: Prentice-Hall, Inc., 1993. 

[37] B. Polyak, Introduction to Optimization. New York: Optimization 
Software, Inc., 1987. 


3.1 


The almost 



